AD  6784 


co 

CD 

Division  of  Engineering  Mechanics 
STANFORD  UNIVERSITY 

Technical  Report  No.  183 
Contract  Report  No.  10 


Division  of  Engineering  Mechanics 


STANFORD  UNIVERSITY 


Technical  Report  No.  183 
Contract  Report  No.  10 


ELASTIC-PLASTIC  DEFORMATION  AT  FINITE  STRAINS 


By 

E.  H.  Lee 


June  1968 
Prepared  under 

Contract  No.  DA-04-200-AMC-659(X) 
Sponsored  through  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground  Maryland 


Thin  Scr:?  ' 
ar.d  eaiu;  u*. 


1  '  r''Io-'3G 

.  -  .i. 


•  a<‘. 


v'  ,  jr-r rtp.'jw-wr 


ELASTIC- PLASTIC  DEFORMATION  AT  FINITE  STRAINS 
E.  H.  Lee 

Stanford  University 


Abstract 

In  some  circumstances,  elastic-plastic  deformation  occurs  in  which 
both  components  of  strain  are  finite.  Such  situations  fall  outside  the 
scope  of  classical  plasticity  theory  which  assumes  either  infinitesimal 
strains  or  plastic-rigid  theory  for  large  strains.  The  present  theory 
modifies  the  kinematics  to  Include  finite  elastic  and  plastic  strain 
components.  For  situations  requiring  this  generalization,  dilatational 
influences  are  usually  significant  including  thermo-mechanical  coupling. 
This  is  introduced  through  the  consideration  of  two  coupled  thermodynamic 
systems:  one  comprising  thermo-elasticity  at  finite  strain  and  the  other 
the  Irreversible  process  of  dissipation  and  absorption  of  plastic  work. 
The  present  paper  generalizes  a  previous  theory  to  permit  arbitrary 
deformation  histories. 

Introduction 

Fig.  1  shows  the  stress-strain  curve  in  tension  or  shear  for  a 
ductile  metal  such  as  aluminum.  Plastic  deformation  sets  in  at  A  ,  and 
the  material  work  hardens  along  ABE  .  Unloading  from  B  is  elastic 
along  BC  ,  and  the  elastic  strain  at  B  ,  CD  in  magnitude,  is  recov¬ 

ered  on  unloading  down  to  zero  stress  at  C  .  The  remaining  strain  OC 
is  the  plastic  strain  associated  with  loading  to  the  point  B  .  The 

magnitude  of  the  elastic  strain  CD  is  equal  to  the  yield  stress  divided 

-3 

by  the  elastic  modulus  and  is  of  the  order  10  .  Plastic  flow  can 
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continue  along  BE  to  strains  of  the  order  unity.  Plasticity  theory 
incorporates  these  features  by  considering  two  limiting  situations 
which  cover  many  practical  circumstances.  If  strains  are  large,  of 
the  order  unity,  as  in  many  metal-forming  problems,  the  elastic  com¬ 
ponent  of  strain  is  considered  negligible  in  comparison  and  plastic- 
rigid  theory  is  utilized.  Since  plasticity  is  an  incremental  or  flow 
type  phenomenon,  strain  rates  based  on  the  velocity  field  appear  in 
the  stress-deformation  relation,  which  with  the  equations  of  motion 
and  boundary  conditions  determines  the  velocity  field  at  each  Instant 
of  time,  and  hence  the  total  deformation  by  integration.  Elastic- 
plastic  theory  has  been  developed  to  consider  situations  in  which  elas¬ 
tic  and  plastic  components  of  strain  are  of  the  same  order,  that  is 
-3 

about  10  as  stated  above,  so  that  infinitesimal  strain  analysis  can 
be  used.  There  are  however  situations  of  plastic  flow  under  high  pres¬ 
sure,  such  as  can  occur  in  explosive  forming,  in  which  the  elastic  and 

plastic  components  of  strain  are  each  finite,  so  that  neither  limit 

* 

covered  in  classical  plasticity  theory  applies  [1]  .  A  theory  incorp¬ 
orating  finite  strain  for  both  elastic  and  plastic  components  is  pre¬ 
sented  below,  which,  in  addition  to  its  appropriateness  for  high 
pressure  loading  problems,  provides  a  uniformly  valid  theory  encompas¬ 
sing  the  two  limits  covered  by  classical  theory. 

As  discussed  in  the  next  section,  the  kinematics  of  finite  strain 
invalidates  the  usual  assumption  that  the  total  strain  is  the  sum  of 
elastic  and  plastic  components.  This  assumption  arises  in  infinitesimal 

Numbers  in  square  brackets  refer  to  the  references  at  the  end  of 
the  paper. 
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straln  analysis  because  of  linearity  and  in  the  one  dimensional  case 
depicted  in  Fig.  1.  Finite  elastic  strains  are  usually  predominently 
dllatational,  since  increase  of  elastic  shear  strain  components  beyond 
the  elastic  limit  (~  10  )  produces  plastic  flow.  The  finite  dilatation 

introduces  significant  thermo-mechanical  coupling  effects  which  must  be 
included  to  provide  a  satisfactory  theory.  With  plastic  flow  this 
demands  the  inclusion  of  irreversibility. 

Kinematics 

As  discussed  in  [2],  the  kinematics  of  elastic-plastic  deformation 

can  be  analysed  on  the  basis  of  the  configurations  depicted  in  Fig.  2. 

Using  a  single  set  of  rectangular  Cartesian  coordinates  (x^,  ,  x^)  , 

particles  in  an  undeformed  body  at  a  uniform  base  temperature  8  have 

o 

positions  represented  by  the  column  vector  x*  .  In  the  deformed  state 

at  time  t  ,  following  elastic-plastic  deformation,  they  occupy  the 
3 

positions  x  ,  according  to  the  mapping 

3  3  1 

x  =  x  (x  ,t)  (1) 

If  each  element  of  the  body  is  unstressed  and  reduced  to  the  base  tem- 

2 

perature,  the  configuration  x  results,  according  to  the  mapping 

2  2  1 

x  -  xZ(x\t)  (2) 

If  a  body  is  reduced  to  the  uniform  base  temperature  0q  and  the  surface 

tractions  are  removed,  a  residual  stress  distribution  will  normally  remain 

following  non-homogeneous  plastic  flow.  To  remove  this  the  body  will 

usually  have  to  be  considered  dissected  into  small  elements  and  the  cor- 

2 

responding  mapping,  (2),  will  be  discontinuous.  The  configuration  x 
at  uniform  base  temperature  and  zero  stress  comprises  the  pure  plastic 
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deformation,  since  thermal  expansion  and  elastic  t,frain  components  are 
both  zero. 

As  discussed  in  [2],  local  deformation  in  the  neighborhood  of  each 
particle  at  time  t  can  be  conveniently  expressed  in  terms  of  the 
deformation  gradient  matrix: 

F  =  axj/axj  (3) 

For  the  case  when  the  mapping  (2)  ii  continuous  and  differentiable,  the 

1  2 

plastic  deformation  from  x  to  x  can  be  expressed  by  the  deformation 
gradient  matrix: 

FP  =  dx^/dXj  (4) 

2  3 

and  the  thermo-elastic  deformation  from  x  to  x  by: 

m  <v 

Fe  =  dx^/dXj  (5) 

The  chain  rule  for  partial  differentiation  then  gives  the  matrix  product 
relation: 

F  =  Fe  FP  (6) 

In  the  case  of  a  discontinuous  mapping  (2),  the  deformation  gradients  (4) 
and  (5)  reduce  to  matrices  of  local  linear  mappings  in  the  limit  of  small 
material  elements,  which  also  satisfy  (6)  but  not  the  continuity  conditions 
for  partial  derivatives.  The  analysis  of  elastic-plastic  deformation  is 
not  essentially  modified  by  this  circumstance. 

The  matrix  product  (6)  provides  the  relation  between  elastic,  plastic 
and  total  deformation  valid  for  finite  strains.  It  replaces  the  usual 
assumption  that  the  total  strain  is  the  sum  of  elastic  and  plastic  compon¬ 


ents  : 
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eij 


(7) 


That  the  latter  is  not  valid  at  finite  strain  is  clear  since  the  plastic 
displacement 


and  elastic  displacement 


are  additive: 


where 


(8) 


The  finite  strain  components  are  non-linear  expressions  in  the  displace¬ 
ments,  and  hence  will  not  be  additive.  Infinitesimal  strain  components 
are  linear  in  the  displacements  so  that  (7)  will  hold  in  that  case.  From 
Fig.  1  it  is  clear  that  (7)  also  holds  in  that  one-component  case. 

Nonlinearity  is  not,  however,  the  most  significant  feature  differ¬ 
entiating  (6)  and  (7).  The  matrix  product  (6)  is  not  in  general  conmuta- 
tive,  which  corresponds  to  the  fact  that  the  plastic  deformation  F*3  is 

a  functional  which  represents  the  past  history  of  plastic  flow  that  has 

£ 

already  occurred,  while  the  elastic  deformation  matrix  F  is  a  function 
of  the  current  stress  and  pre-multiplies  F*5  .  For  finite  strain  the 
elastic  and  plastic  deformations  are  coupled  through  (6)  in  a  much  more 
involved  manner  than  that  corresponding  to  the  addition  (7). 

3  2 

The  un-stressing  process  relating  the  configurations  x  and  x 

is  not  uniquely  defined  since  an  arbitrary  rigid  body  rotation  can  be 
2 

superposed  on  x  and  still  leave  that  configuration  unstressed.  Since 
2 

in  general  the  x  configuration  is  not  continuous,  such  an  arbitrary 
rotation  can  differ  for  each  element  of  the  body. 


m  - 
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Since  we  shall  assume  that  the  constitutive  relation  for  elastic- 
plastic  deformation  contains  the  deformation  history  through  and 

G 

F  only,  and  does  not,  for  example,  include  strain  gradient  or  size 
effects,  we  need  consider  only  homogeneous  states  of  stress  and  defor¬ 
mation  in  formulating  the  constitutive  law.  Unloading  the  surface  of 

a  body  so  deformed  would  then  reduce  the  stresses  to  zero  at  uniform 

2 

temperature,  and  no  questions  of  a  discontinuous  configuration  x 
arise. 

Work  Hardening 

Plastic  flow  is  an  extremely  complicated  phenomenon,  and  a 
theory  which  attempts  to  include  all  its  aspects  is  likely  to  be  pro¬ 
hibitively  difficult  to  either  formulate  or  apply.  We  shall  therefore 
develop  an  approach  based  on  the  classical  theory  of  Isotropic  work 
hardening  generalized  to  include  additional  features  which  finite  strain 
analysis  appears  to  require.  Other  features  not  considered  in  this 
approach,  such  as  a  BauBchinger  effect  or  the  localized  modification  of 
the  yield  surface  due  to  plastic  flow,  could  be  incorporated  as  they  are 
in  developments  of  classical  plasticity  theory. 

Our  analysis  will  follow  the  approach  to  isotropic  work  hardening 
ably  expounded  by  Hill  ([3]  Chap.  II)  for  strain  increments  with  linear 
elastic  law.  The  system  of  yield  surfaces  associated  with  continued 
work  hardening  consists  of  a  one  parameter  family  of  non-intersecting 
surfaces  in  stress  space.  Plastic  flow  generating  a  stress  corresponding 
to  a  point  on  one  of  these  surfaces,  determines  that  surface  as  the  cur¬ 
rent  yield  surface  until  further  plastic  flow  takes  place.  This  yield 
surface  is  independent  of  the  history  of  loading  as  long  as  the  stress 
has  not  previously  exceeded  this  yield  value.  The  yield  condition  is 
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found  to  be  effectively  independent  of  superposed  hydrostatic  pressure, 
as  is  the  hardening  expressed  by  the  stress-strain  relation  when  stress 
deviator  components  are  plotted.  The  yield  limit  as  plastic  flow  develops 
is  found  to  be  a  function  of  the  expenditure  of  work  associated  with  the 
plastic  strain.  Since  superposed  hydrostatic  pressure  is  found  not  to 
influence  the  yield  limit  or  the  work  hardening  phenomenon,  it  follows 
that  no  volume  change  occurs  in  plastic  deformation,  for  otherwise  the 
superposed  hydrostatic  pressure  would  contribute  to  the  work  of  plastic 
flow  and  hence  influence  plasticity.  This  deduction  is  found  to  be  in 
good  agreement  with  experiment.  At  finite  elastic  strains,  thermo- 
meciianical  coupling  influences  become  significant,  so  that  temperature 
variation  must  be  considered  in  the  theory,  and  thermo-elastic  theory 
utilized  for  the  elastic  component  of  deformation.  For  consistency, 
thermo-mechanical  coupling  Influences  must  also  be  Incorporated  in  the 
plasticity  relations. 

To  formulate  a  theory  on  the  basis  of  the  features  mentioned  above, 
one  must  separate  che  work  stored  in  elastic  strain  energy  from  the  total 
work  expended  in  plastic  flow.  For  the  infinitesimal  strain  case  such 
separation  follows  at  once  using  (7).  The  rate  of  expenditure  of  work 
on  the  material  per  unit  volume  is  given  by  : 


clearly 


aij 


+ 


w*  -  .  ilj 


(9) 
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gives  the  rate  of  expenditure  of  work  stored  in  elastic  deformation, 
while 

•  P  .p 

w 

gives  the  rate  of  increase  of  plastic  work  which  governs  the  work 
hardening  phenomenon.  An  analogous  separation  must  be  achieved  when  (6) 
replaces  (7)  for  the  finite  deformation  case.  This  was  effected  in  [2] 
only  for  restricted  deformations  in  which  principal  directions  remained 
fixed  in  the  body  throughout  the  deforming  process. 

3  1 

Consider  a  body  deforming  through  the  configuration  x  (x  , t)  shown 
in  Fig.  2.  The  rate  of  work  expended  in  deforming  the  material  of  the 
body,  that  is  the  total  rate  of  work  less  the  rate  of  generation  of 
kinetic  energy,  is  given  by  (see,  for  example  [4]  p.  87) 


W 


9vi  3 
aii  “J  dV 
V3  j  H 


(10) 


where  v^  is  the  particle  velocity,  and  V  the  volume  in  the  configura¬ 


tion  x 


and 


It  will  prove  convenient  to  use  matrix  notation  and  write 


T  =  a 


3 

3 

„  1 

9vi 

i 

cv 

< 

!-*• 

9xk 

,  3 

„  1 

3 

SXj 

3xk 

ij 


=  F  F  =  L 


(ID 


(12) 


F  being  defined  in  (3).  F  represents  time  differentiation  at  constant 
x^  ,  and  hence  convected  differentiation  at  a  particular  particle,  since 
the  body  is  at  rest  in  the  initial  configuration  x^  .  In  order  to 


mo  i  jfa. 
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determine  the  work  expended  on  a  specified  element  of  material,  it  is 
convenient  to  change  the  variables  to  integration  over  the  body 
configuration  in  .  Equation  (10)  then  becomes: 


W  =  f  tr(T  L)  det(F)  dV1 

j 


(13) 


where  "tr"  stands  for  the  trace  of  the  matrix  argument  (the  sum  of 
diagonal  terms)  and  "det"  for  determinant.  Det(F)  is  the  Jacobean  of 
the  transformation  of  variables  for  the  volume  integral.  Substituting 

(6)  for  F  in  (13)  and  (12),  and  using  the  product  formula  for  differen¬ 
tiation  of  a  matrix  product,  (13)  becomes: 

-1  -l 

W  =  [  tr[T(F®FP  +F®FP)FP  F®  ]  det(F®)det (Fp)  dV1  (14) 

J  ..1 . 


Zero  volume  change  in  plastic  flow  determines: 

det(FP)  =  1 


(15) 


and  hence  W  divides  into  the  two  components  which  will  be  tentatively 
labeled  elastic  and  plastic  work  rates: 

■1 


W®  -  [  tr[TF®  F®  ]  det(FC)  dV1 

j  ~~  •- 


(16) 


and 


-1  -1 


WP  =  J  tr[TF®  FP  FP  F®  ]  det(F®)  dV1  (17) 


According  to  the  kinematics  discussed  in  connection  with  Fig.  2, 

6  1 

F  (x  ,t)  represents  the  deformation  gradient  associated  with  elastic  strain 

of  the  plastically  deformed  body  at  each  instant.  Thus  (16)  has  the  form  of 

the  rate  of  storage  of  elastic  strain  energy  in  an  elastic  body  taken  through 

0 

the  deformation  history  F  (x  ,t)  .  However,  care  must  be  exercised  in 
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transferring  this  interpretation  to  the  elastic-plastic  case,  since  the 

2 

unstressed  state  x  is  continuously  deforming  plastically,  whereas  norm¬ 
ally  for  purely  elastic  deformation  the  unstressed  state  is  a  fixed  undls- 
torting  configuration  of  the  body.  Since  the  elastic  strain  energy  depends 
on  the  elastic  deformation  gradient  and  the  elastic  characteristics  of  the 
material,  (16)  might  not  represent  the  rate  of  increase  of  elastic  strain 
energy  if  the  elastic  characteristics  are  changing  due  to  the  continuing 

plastic  flow.  For  example,  the  elastic  strain  energy  could  change  at  con- 
e 

stant  F  in  that  case.  Thus  the  interpretation  of  (16)  and  (17)  depends 
crucially  on  the  coupling  between  elastic  and  plastic  deformations. 

We  have  already  seen  that  characteristics  of  plastic  flow  are  insen¬ 
sitive  to  hydrostatic  pressure  and  so  to  the  elastic  dilatation.  Corres¬ 
pondingly,  elastic  characteristics  are  insensitive  to  plastic  flow.  The 
usual  influence  assumed  of  plastic  flow  is  on  the  yield  limit  with  the 
elastic  constants  remaining  unchanged.  Hill  ([3]  p.  24)  states  that  appre¬ 
ciable  anisotropy  only  appears  at  large  strains.  Indeed,  stress  analysis 
is  usually  based  on  standard  elastic  constants  even  though  the  components 
may  have  been  cold-worked.  Moreover,  to  our  knowledge,  little  information 
exists  on  the  modification  of  elastic  constants  due  to  plastic  flow.  This 
state  of  affairs  is  consistent  with  our  knowledge  of  the  physical  basis  of 
plasticity  associated  with  the  generation  and  migration  of  dislocations. 
Even  in  a  highly  plastically  deformed  material  the  basic  crystal  structure 
is  retained  with  disturbance  from  the  regular  lattice  of  only  a  very  small 
proportion  of  the  atoms.  Thus  the  elastic  constants  and  heat  capacity  are 
not  appreciably  influenced  by  plastic  flow.  We  will  therefore  assume  the 
elastic  characteristics  to  be  governed  by  an  invariant  isotropic  thermo¬ 


elastic  law. 
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3 

Thus  the  stress  X  In  the  configuration  x  will  depend  only  on 
e 

F  and  temperature  9  ,  since  plastic  deformation  and  rigid-body  rotation 

2 

of  the  unstressed  state  x  ,  will  not  modify  the  stress  for  an  isotropic 
elastic  material  with  invariant  elastic  properties.  The  rate  of  investment 
of  elastic  energy  in  the  thermo-elastic  component  of  the  system  will  thus 
be  the  same  as  that  for  a  purely  elastic  body  deformed  by  Fe  at  tempera¬ 
ture  0  from  a  fixed  undistorted  state.  Equation  (16)  thus  represents 
“.he  rate  of  expenditure  of  thermo-elastic  energy  invested  in  the  body 

mechanically.  The  elastic  properties  can  be  conveniently  expressed  in  terms 

T 

6  6  6  C 

of  the  Helmholtz  free  energy  per  unit  mass  i|f(C  ,0)  ,  where  C  =  F  F 
(superscript  T  denotes  transpose  of  the  matrix) ,  and  the  stress  and  entropy 
in  the  thermo-elastic  system  are  then  given  by  ([5]  p.  302  and  309,  or  [6] 
p.  199): 

T  =  2p  Fe  F®  /det(Fe>  (18) 

-  0  '  3Cc  ' 

se  =  .  ai 

ae  (19) 

where  n  is  the  initial  density,  and  Se  the  entropy.  Examples  of  forms 
r  0 

of  the  free  energy  function  f  are  given  in  [2]  and  [7]. 

It  is  interesting  to  note  that  the  lack  of  coupling  between  the  laws 

of  elastic  and  plastic  deformation  described  above  is  also  manifest  in 

the  incompressible  property  of  plastic  flow.  For  a  change  in  density  of 

the  unstressed  state  would  imply  modified  elastic  characteristics. 

Since  (16)  gives  the  rate  of  elastic  working,  (17)  must  represent  the 

mechanical  work  being  expended  in  plastic  flow.  This  appears  to  exhibit 

£ 

involved  coupling  with  the  elastic  deformation  F  ,  but  this  can  be  largely 
eliminated  as  shown  below.  Since  we  shall  be  concerned  with  stress-strain 
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relations,  as  mentioned  earlier,  we  can  without  loss  assume  homogeneous 
stress  and  deformation  distributions  and  replace  the  volume  integral  in 
(17)  by  wP  ,  the  rate  of  plastic  work  per  unit  undeformed  volume. 


wP  =  tr(TFe  FP  FP  Fe  )  det(F6)  (20) 

2 

It  was  pointed  out  in  the  discussion  on  kinematics  that  the  state  x 

is  not  uniquely  defined  since  de-stressing  can  involve  arbitrary  rotation. 

-1 

Analogously  with  (12),  the  terra  FP  Fp  represents  the  velocity  gradient 

2 

of  plastic  flow  in  x  ,  and  its  symmetric  part  gives  the  velocity  strain, 

2  2 

often  loosely  called  the  rate  of  strain,  in  x  .  When  x  is  arbitrarily 
rotated,  the  work  associated  with  this  rate  of  plastic  strain  can  clearly 
be  completely  changed.  For  example,  stretch  in  the  direction  of  a  tensile 
stress  could  become  lateral  stretch  by  a  90°  rotation  of  x^  ,  with  conse¬ 
quent  drastic  modification  of  the  plastic  work  term.  Since  we  assume  the 
material  to  be  isotropic,  the  plastic  strain  rate  tensor  will  he  parallel 

to  the  stress  tensor,  and  the  work  relation  as  it  actually  takes  place  will 

2 

be  achieved  by  considering  x  to  be  defined  by  de-stressing  without  rota- 

tion.  This  special  choice  of  F  ,  denoted  by  F  ,  will  then  be  symmetric 

and  represent  pure  stretch  in  the  directions  of  the  principal  axes  of  stress, 

2 

in  conformity  with  the  isotropic  elastic  assumption.  The  state  x  is 

Q 

then  uniquely  prescribed.  The  availability  of  this  choice  for  F  can  be 

seen  analytically  since  the  polar  decomposition  theorem  ([5]  p.52)  permits 

e 

any  matrix  of  the  form  F  to  be  expressed  as: 


F6  = 


e  e 
V  R 


(21) 


w 
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Q  6 

where  V  is  symmetric  and  R  a  rotation  matrix.  Thus  (6)  can  be 
written 


ii 


F  -  FeFP  =  VeReFP  =  Ve(ReFP)  =  FeFP 


(22) 


Ve  and  FP  =  ReFp 


where  F 

Since  the  trace  of  a  product  of  matrices  is  not  influenced  by  cyclic 

permutation  of  the  elements  as  is  clear  from  the  subscript  representation, 

—2 

(20)  can  be  written  for  the  special  choice  x  : 

-1  .  -1 

w15  =  tr(Fe  T^FPp  )  dettf6)  (23) 

The  trace  is  a  scalar  invariant  of  the  tensor  arguments,  so  that  change  of 

axes  does  not  modify  its  value.  Selecting  axes  parallel  to  the  directions 

e 

of  principal  stress,  T  and  F  become  diagonal  matrices  which  are  commu- 
tative,  so  that  F  and  F  cancel,  and  (23)  reduce-  fo 


-1 

wP  =  tr (T  F  P  )  det(p) 


(24) 


which  retains  this  form  on  a  transformation  back  to  the  original  coordinates. 
This  can  alternatively  be  written 


^—2 

^  =  tI  det(Fe) 


(25) 


—2  —2 
where  v^  is  the  velocity  of  plastic  straining  in  the  configuration  x  . 

Thus  the  plastic  rate  of  work  term  takes  on  the  simple  form  of  work  done 

by  the  stress  T  on  the  plastic  rate  of  deformation  with  a  scalar  coupling 

term  to  the  elastic  deformation  :  det(Fe)  .  It  is  now  seen  that  the 

result  is  identical  with  that  discussed  in  [2]  in  which  the  analysis  was 


1rTJHiMMinrr  ir*  — — 
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limited  to  total  deformation, throughout  the  loading,  having  principal 

directions  fixed  in  the  body,  for  wSich  all  matrices  considered  were 

thus  diagonal,  so  that  the  matrix  contribution  of  F  cancelled  directly. 

In  retrospect,  since  the  rate  of  total  work  (13)  and  the  rate  of 

2 

elastic  work  (16)  do  not  depend  on  the  particular  choice  of  x  ,  nor  can 
the  rate  of  plastic  work.  Substituting  (21)  for  Fe,  and  permuting  the 
elements  in  the  trace,  (17)  becomes 

WP  -  J  tr(V6  TVCReFP  FP  R6  )  det  (F®)  dV1  (26) 

V1 

e 

For  isotropic  elastic  behavior  the  rotation  R  prior  to  the  pure 
£ 

stretch  V  in  (21)  does  not  affect  the  stress  T  which  is  parallel  to 

e  e  ^  e 

the  pure  stretch  tensor  V  .  Thus  the  V  and  V  terms  cancel. 

-1  ~ 

e  e 

The  R  ,  R  terms  simply  represent  rotation  of  the  plastic  velocity 

gradient  tensor,  or  equivalently  velocity  strain  since  T  is  symmetric. 

2 

Thus,  (26)  tells  us  that,  for  arbitrary  x  ,  we  must  rotate  the  plastic 
velocity  strain  term  by  the  elastic  rotation  to  determine  the  plastic 

-2 

rate  of  work,  and  this  is  equivalent  to  making  the  special  selection  x 
Having  determined  the  rate  of  plastic  work,  we  can  formulate  an 
isotropic  work  hardening  law  analogous  to  that  described  by  Hill  ([3] 

Chap.  II).  However,  since  thermal  effects  can  be  significant  in  the  present 
study,  the  influence  of  varying  temperature  on  the  yield  limit  should  be 
included,  and  we  thus  need  to  generalize  Prager's  analysis  for  the  rigid- 
plastic  case  [8].  This  was  done  in  [2]  for  the  special  case  of  principal 
directions  fixed  in  the  body. 

The  rate  of  plastic  work  invested  per  unit  initial  volume,  (24),  can 
be  written: 

-1 

wP  =  tr[(T  det  T)  F9?9  ] 


(27) 
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and  thus  T  can  be  replaced  by  its  deviator  without  modifying  (27)  or 
(28).  Thus  hydrostatic  pressure  only  influences  the  plastic  work  through 

g 

the  scalar  factor  det(F  )  ,  which  is  a  term  expressing  geometrical  non¬ 
linearity  associated  with  reduction  in  scale  of  the  initially  unit  volume 
element  and  hence  also  of  the  areas  over  which  the  stress  acts  and  the 
velocity  through  which  the  forces  are  working.  As  mentioned  above,  the 
work-hardening  produced  by  plastic  flow  is  not  appreciably  modified  by  the 
existence  of  superposed  hydrostatic  pressure,  and  physically  it  is  associ¬ 
ated  with  the  generation  of  Increased  dislocation  density.  The  production 
of  a  prescribed  dislocation  density  associated  with  a  prescribed  state  of 
hardening  is  likely  to  require  more  plastic  work  under  high  pressure  rather 
than  less,  since  the  small  coupling  existing  relates  to  a  slight  increase 
of  volume  due  to  disturbance  of  the  regular  crystal  lattice.  To  avoid  the 
contradiction  of  less  work  according  to  (27)  under  hydrostatic  pressure, 
we  suggest  that  the  yield  condition  for  finite  elastic  strain  should  take 
the  form 

f[T  (det  ?*)]  =  c  (31) 


0MH&X***  — 
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where  f  takes  on  one  of  the  usual  forms,  such  as  Mises  or  Tresca,  but  could 

be  a  more  general  function  of  ^  and  ,  the  second  and  third  invariants 

of  the  argument  deviator,  as  discussed  in  [3].  Note  that  det(Fe)  is  a 

scalar  multiplier  which  does  not  alter  the  orientation  of  the  tensor  argument. 

The  yield  condition  (31)  reduces  to  classical  yield  conditions  in  the  case  of 

infinitesimal  elastic  strain  for  which  det(F  )  can  be  taken  equal  to  unity. 

The  static  experiments  carried  out  to  date  from  which  the  form  of  (31)  might 

have  been  tested,  are  not  sufficiently  precise  to  achieve  this  in  view  of 

£ 

the  closeness  of  det(F  )  to  unity. 

For  a  work  hardening  material  subject  to  temperature  changes,  c  in 
(31)  has  been  considered  to  be  a  function  of  w^  and  the  temperature  0 
[2].  However,  it  is  known  that  increasing  the  temperature  increases  the 
mobility  of  dislocations  and  hence  permits  the  generation  of  increased  dis¬ 
location  density,  and  so  of  hardening,  with  a  smaller  expenditure  of  plas¬ 
tic  work.  Hence  a  functional  rather  than  a  function  of  w*5  is  needed  in 
the  expression  for  c  in  (31),  which  represents  the  hardening,  and  the 
form 


c 


4 


t 

aid)  dT»9] 


(32) 


suggests  itself.  The  term  0(9)  will  be  an  increasing  function  of  0  to 
achieve  the  desired  temperature  -  dislocation  mobility  influence.  Relation 
(31)  with  changing  c  according  to  (32)  will  yield  the  desired  set  of 
similar  yield  surfaces  when  f  is  an  isotropic  homogeneous  function  of 
its  argument  components,  as,  for  example,  is  the  Mises  condition. 

Following  Prager  [8]  it  is  revealing  to  write  the  yield  condition 


(31),  (32)  in  the  form: 
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where  Q  stands  for 


f (Q)  -  c(cp,0)  <  0 


9  =  T(det  F*) 


and  cp  represents  the  functional 


J  o(9)  wF 


The  equality  sign  in  (33)  denotes  plastic  flow  or  neutral  loading,  and  the 
inequality  sign,  an  elastic  state.  For  continuing  plastic  flow,  differen¬ 
tiation  of  (33)  with  the  equality  sign  gives 


&lQ-a£.i.&£g.o 

99  9cp  ™  99 


The  positive  nature  of  plastic  work  and  its  hardening  effect  determines 

£*>0 

so  that  continued  plastic  flow  is  associated  with 

If  ?  -  If e  > 0  (37) 

and  this  becomes  the  loading  condition.  Note  that  in  the  isothermal  case 

•  » 

this  simply  corresponds  to  the  usual  condition  f  >  0  . 

In  order  to  define  the  law  governing  rate  of  plastic  flow,  it  is  con¬ 
venient  to  work  in  terms  of  the  plastic  velocity  strain,  often  loosely  called 
the  rate  of  strain: 


2  -2 

9v,  dv.N 


D9  =  (I?  +  L9  )/2  -  +  -^0/2  ;  L9  =  Ff 


Because  T  is  symmetric,  (27)  can  be  written: 

wP  -  tr[(T  det  F*)  D9] 


mz- w«suw»**»"fc 
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Whereas  this  relation  combined  with  (35)  and  (36)  was  sufficient  to 
determine  in  the  case  of  one  dimensional  strain  [2],  since  only  one 

non-zero  component  remained  having  introduced  incompressibility,  in  gen¬ 
eral  we  need  the  concept  oi  plastic  potential  to  determine  the  several 
components.  This  takes  the  form  [8] 

5P  =  ktf(?f§  «>  (40) 

where  k  is  a  function  of  the  current  state  to  be  determined.  The  term  in 
the  parentheses  gives  the  rate  of  loading  in  analogy  with  (37),  and  the 
normality  condition  of  the  strain  rate  components  and  the  yield  surface 
according  to  the  tensor  term  Qf/^Q  ,  leads  to  a  maximum  rate  of  work 
principal  and  a  uniqueness  of  solution  argument  as  for  classical  plasticity. 

Since  the  equality  sign  applies  in  (33)  we  can  substitute  (40)  into 
(36),  using  (39)  and  obtain: 


If  ^  -  If 4  -  U“<0)  "M  k  H  ^ 


(41) 


The  scalar  quantities  k  and  (^  Q  -  9)  factorize  from  the  trace 

expression,  the  latter  simplifies  by  use  of  Euler's  theorem  since  f  is 

homogeneous  in  Q  ,  and  (41)  reduces  to 

1  =  k  a  n  f  (42) 

dcp 

where  n  is  the  order  of  Q  in  f(Q)  ,  and  (40)  becomes 


?  ■  If  <H  s  -  If  'e)/a  n  £  I  <43) 

This  relation  is  equivalent  to  that  for  isothermal  work  hardening  given  by 
Hill  ([3]  p.  38)  and  is  seen  not  to  be  appreciably  complicated  by  the  inclu¬ 
sion  of  finite  strain,  non-linear  elasticity  and  thermo-mechanical  coupling. 


The  Thermodynamic  System 


The  analysis  presented  already,  permits  the  determination  of  deforma¬ 
tion  given  the  stress  history  apart  from  the  temperature  variation.  To 
determine  this, thermo-mechanical  coupling  in  plastic  flow  must  be  consid¬ 
ered  in  addition  to  that  included  in  the  thermo-elastic  theory  (18)  and 
(19).  The  approach  follows  directly  that  used  in  the  fixed  principal 
directions  case  [2].  We  have  seen  that  the  elastic  and  plastic  phenomena 
are  only  loosely  coupled,  and  that  thermo-elastic  properties  depend  on  the 
basic  crystal  structure  whereas  plastic  flow  depends  on  generation  and 
migration  of  dislocations  which  leave  most  of  the  crystal  structure  intact. 
We  thus  base  the  thermodynamics  on  two  interacting  systems  i)  reversible 
thermoelasticity  governed  by  (18)  and  (19)  and  the  corresponding  specific 
heat  relations  and  ii)  irreversible  plastic  flow,  comprising  some  dissi¬ 
pation  of  the  plastic  work  to  heat  which  appears  as  a  source  in  the  thermo¬ 
elastic  system  and  hence  as  an  increase  in  the  elastic  entropy  in  adiabatic 
loading.  The  randomness  of  the  dislocation  distribution  also  leads  to  an 
increase  in  entropy  associated  with  plastic  distortion.  A  strain  energy 
field  surrounds  each  dislocation  which  absorbs  some  of  the  work  of  plastic 
flow  as  the  dislocation  density  grows.  For  the  specific  thermo-mechanical 
coupling  process  generated  by  plastic  flow,  we  adopt  that  measured  by 
Taylor  and  Farren  [9]  and  Taylor  and  Quinney  [10],  in  which  temperature 
rise  was  observed  in  elastic-plastic  tests  which  exceeded  that  expected 
from  thermo -elastic  theory  by  an  amount  equivalent  to  a  heat  source  equal 
to  ywP  per  unit  initial  volume,  where  y  is  a  factor  which  varies  slowly 
from  cbout  0.9  to  unity  with  increasing  plastic  flow.  The  rest  of  the 
plastic  work  rate,  (l-y)wP  ,  is  considered  to  be  stored  in  the  developing 
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dislocation  system.  Thus  the  internal  energy  per  unit  mass  associated  with 
thermo-elastic  deformation,  E  ,  and  that  stored  in  the  dislocation  dis¬ 
tribution  Ep  ,  are,  for  adiabatic  loading,  determined  by 

p^E  =  w  +  yw  (44) 

and 

PqEp  =  (1-y)  wP  (45) 

0 

Growth  of  elastic  entropy  S  arises  due  to  the  mechanical  energy  dissi¬ 
pated  in  plastic  flow,  and  is  given  by 


We  shall  consider  y  to  take  on  the  constant  value  0.9.  These  relations 
assume  continuous  variations  of  deformation.  In  the  case  of  shock-wave 
loading,  involving  discontinuous  changes  of  stress  and  strain,  additional 
increase  in  entropy  occurs  over  and  above  that  given  by  (46). 

Combination  of  the  contributions  developed  in  this  and  the  previous 
section  determines  the  total  stress-deformation  relation.  Equation  (18) 
for  elastic  strain,  (33)  and  (43)  for  plastic  flow  and  (19),  (44)  and  (46) 
for  the  coupling  between  these  components. 

Discussion 

The  rather  concise  structure  of  the  theory  developed  in  this  paper 
arises  because  of  the  comparatively  loose  coupling  between  the  elastic  and 
plastic  phenomena.  The  coupling  needed  to  achieve  simplicity  and  the  reten¬ 
tion  of  the  basic  structure  of  classical  theory  was  determined  by  the 
expression  for  plastic  work  rate  (27)  which  led  to  the  appropriate  yield 
relation  (31).  The  introduction  of  the  influence  of  plastic  flow  on  elastic 
properties  would  clearly  lead  to  a  more  elaborate  theory,  but  fortunately 
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these  effects  are  not  dominant. 

With  the  present  development  In  mind,  it  is  Interesting  to  review  the 
presentation  by  Hill  (  [2]  Chap.  II).  He  carefully  limits  (7)  to  incre¬ 
ments  since  he  is  concerned  with  finite  plastic  deformation  and  infini¬ 
tesimal  elastic  strains.  His  strain  increment  relation  corresponds  to 
our  velocity  strain  relation  from  (12) 

D  =  (L  +  LT)/2  (47) 

where 

L  =  FeF6  +  FeFP  FP  F6  (48) 

Clearly  the  elastic  and  plastic  velocity  strains  components 

T.e  =  FeFe  ,  LP  =  frPFP  (49) 

are  not  additive  to  give  the  total  velocity  strain, but  for  infinitesimal 
£ 

elastic  strains  F  ~  I  ,  the  unit  matrix,  and  additivity  applies  to  some 
order  of  approximation.  For  the  fixed  principal  directions  case,  the 
matrices  in  (48)  are  diagonal  and  the  Fe  terms  cancel,  leading  to  additi¬ 
vity,  as  utilized  in  [2]. 
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Fig.  1.  Stress-strain  curve  for  a  ductile  metal 


Fig.  2.  The  kinematics  of  elastic-plastic  deformation 
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